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Temperature and magnetic field studies of the elastic constants of the chromium spinel CdCr2 04 
show pronounced anomalies related to strong spin-phonon coupling in this frustrated antiferromag- 
net. A detailed comparison of the longitudinal acoustic mode propagating along the [111] direction 
with a theory based on an exchange-striction mechanism leads to an estimate of the strength of the 
magnetoelastic interaction. The derived spin-phonon coupling constant is in good agreement with 
previous determinations based on infrared absorption. Further insight is gained from intermediate 
and high magnetic field experiments in the field regime of the magnetization plateau. The role of 
the antisymmetric Dzyaloshinskii-Moriya interaction is discussed. 

PACS numbers: 43.35.-(-d, 72.55.-l-s 



I. INTRODUCTION 

Highly frustrated antiferromagnets (AFs) are of great 
interest due to their potential to realize various uncon- 
ventional phases, even in the classical limit. An example 
in three dimensions is the much studied pyrochlore lat- 
tice, which consists of corner-sharing tctrahedra. On the 
theoretical side, it was noticed quite early that the clas- 
sical Heisenberg AF with nearest-neighbour interactions 
on a pyrochlore lattice has a very unconventional ground 
state, often referred to as a classical spin liquid or a co- 
operative magnet^'^ that is characterized by power-law 
(dipolar) spin correlations.'^'^ Various perturbations in 
such a highly correlated paramagnet are known to have 
a drastic effect on the ground state. The discovery of 
several materials of this kind has led to a strong interest 
in magnets with a pyrochlore lattice. 

A classic example of a pyrochlore Heisenberg AF is 
the spinel compound CdCr204. Here, the Cr^+ ions 
with spin S = 3/2 form the pyrochlore lattice. Al- 
though the Curie- Weiss temperature is Qcw ~ — 70 K 
in this compound, the AF ordering sets in only at Tn ~ 
7.8 K, indicating a high level of magnetic frustration.^ 
The intermediate regime, also referred to as a coopera- 
tive paramagnet, already exhibits strong but short-range 
spin correlations.^ The magnetic ordering at Tn is ac- 
companied by a structural change with loss of inver- 
sion symmetry from the cubic (FdSm) to the tetrag- 
onal (lAi/amd) (Ref. 7) structure. This first order 
magneto-structural phase transformation is caused by a 
spin Jahn-Teller effect, whereby the spins relieve their 
frustration by distorting the crystal.® Neutron-scattering 
studies^'^ show a long-pitched incommensurate coplanar 
spiral spin configuration below Tjv [ordering wave vector 
Q = 27r(0, S, 1), S = 0.09]. This, in turn, is attributed to 



weak Dzyaloshinskii-Moriya (DM) interactions"'^" that are 
present in this compound and assert themselves in the 
magnetically ordered phase. Electron Spin Resonance 
(ESR) measurements suggest^^ that, in a magnetic field 
{H w 5.7 T), the spiral transforms to a four-sublattice 
canted structure. Pulsed-field magnetization data ex- 
hibit a broad plateau between 28 and about 60 T,^'^'^ 
which corresponds to 1/2 of the saturation magnetiza- 
tion of the Cr'^"'" magnetic moments. The fully polar- 
ized state is achieved above 90 T.^-^ The anisotropy of 
the magnetic exchange interactions otherwise appears to 
be negligibly small in this compound.^ All experiments 
clearly indicate the importance of spin-lattice coupling 
in this compound and, indeed, in the whole family of 
chromium spinels (ACr204, A = Zn, Cd, Hg). 

Ultrasound investigation is a powerful method to 
probe spin-lattice interactions and possible lattice 
instabilities-'^. Indeed, it is an ideal experimental tool 
to investigate the chromium spinels as the magnetodis- 
tortive transition is precisely about the interplay of spin 
and lattice degrees of freedom. We have performed mea- 
surements of the relative change of the sound velocity 
in CdCr204 using a phase-sensitive detection technique, 
based on a standard pulse-echo method with a setup as 
described in more detail in Refs. [14,15]. This technique 
is available as well for non-destructive pulsed magnetic 
fields, extending the parameter range for measuring the 
sound velocity and sound attenuation up to very high 
magnetic fields. The measurement accuracy for a rela- 
tive change of sound velocity is of the order of 10~^ for 
static-field measurements and 10^^ for the pulsed-field 
experiments. The sound velocity v is related to the elas- 
tic constant c = pv^ , where p is the mass density of the 
crystal. In this paper, we report on the relative change of 
the sound velocity for a longitudinal acoustic, cl, mode 
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propagating along the [111] direction of CdCr204 (of the 
undistorted cubic crystal) both in static fields up to 18 T 
and pulsed fields up to 63 T. Both magnets were equipped 
with ^He-fiow cryostats. Wide frequency-range piezoelec- 
tric film transducers glued on the parallel surfaces of the 
sample with a two-component epoxy were used for the 
excitation and detection of the longitudinal sound waves. 

The acoustic mode under consideration corresponds to 
cl = l/3(cii + 2ci2 + 4C44), with being the elastic 
constants for a cubic crystal (we have used the familiar 
Voigt notations). The magnetic field was applied along 
the [111] direction, i.e., parallel to the wave vector k and 
the polarization u of the ultrasound wave. The magne- 
tization has been measured using a commercial SQUID 
(superconducting quantum interference device) magne- 
tometer equipped with a 7 T magnet. A high-quality 
single crystal of CdCr204 was obtained by spontaneous 
crystallization from Bi203-V205 flux. The sample thick- 
ness along the [111] direction is 0.83 mm. 

The sound velocity shows anomalous dependences as 
a function of both temperature and magnetic field. We 
provide a theoretical framework to understand these ef- 
fects on the basis of spin-lattice coupling. Our paper 
reveals a fascinating interplay of spin and lattice degrees 
of freedom in this highly frustrated magnet. 

The paper is organized as follows. In Sec. II, we in- 
troduce the microscopic Hamiltonian that captures the 
spin-phonon interactions in CdCr204 and this forms the 
basis of our analysis. The effect of the spin-lattice in- 
teraction on the sound velocity is divided up according 
to the magnetic phases. We start by discussing the be- 
havior in the paramagnetic phase in Sec. III. Then, we 
briefly, discuss the features of the magnetoelastic tran- 
sition in Sec. IV. This is followed by the discussion of 
the low-temperature phase in Sec. V and the effect of 
the magnetic field in Sec. VI. Finally, we summarize our 
conclusions in Sec. VII. The details of various calcula- 
tions are given in Appendices A and B. 



II. THE SPIN-PHONON HAMILTONIAN 

In CdCr204, the CrOe octahedra build an edge- 
sharing network with Cr^+ ions forming a pyrochlore 
lattice. The octahedral symmetry of the crystal field 
splits the five Cr d-orbitals and lowers the energy of the 
three t2g orbitals compared to the two e^ orbitals. Strong 
Hund's coupling aligns the three electrons in the t2g or- 
bital leaving a net spin of 5 = 3/2 on each Cr'^"'' ion. 
The orbital part of the electron wave function forms a 
singlet, and the orbital degrees of freedom are effectively 
quenched. At low temperatures, a spin-only Heisenberg 
Hamiltonian suitably describes the system. The non- 
magnetic Cd ions control the Cr-Cr distance and thereby 
the value of the exchange strength. Further-neighbour 
exchanges are weak owing to the arrangements of differ- 
ent relevant overlapping orbitals. ^^'^^ Thus, the principal 



part of the spin Hamiltonian given by 

i/,'p = ^J.,S..S,-^h.S, (1) 

encodes the antiferromagnetic {Jij > 0) exchange be- 
tween nearest neighbours (ij). The second term denotes 
the usual Zceman coupling to an external magnetic field 
H, with h = 9(UbH {g and are the gyromagnetic ratio 
and the Bohr magneton, respectively). 

In the absence of orbital degrees of freedom, the mag- 
netoelastic coupling is mediated by the dependence of Jij 
on the position of the Cr spins. (Throughout our calcu- 
lations, we treat the spins as classical.) By expanding 
Jij in Eq. (1) around the equilibrium positions of the Cr 
ions to harmonic order, we get 

dJij 1 9^Jr 

■Jij ^ Ja+ ■ Rij + ■^'Rij ■ ■ (2) 

Jo is the equilibrium exchange coupling and Ry = 
Hi — Rj , where is the displacement of i*'* Cr ion from 
its equilibrium position. For the full spin-phonon Hamil- 
tonian, we must introduce phonons. Writing bosonic 
creation and annihilation operators, 6k and b]^, for the 
phonons, we obtain 

H = Hph + Hsp + iJph_sp, (3) 

with 

JJph = hJ2^'^bibi,, (4a) 

k 

Hsp = Jo^S, -Sj-^h-S,, (4b) 

Fph-sp = H1+H2. (4c) 

We are interested in the longitudinal phonon mode with 
u!^ and ek as the bare frequency and the polarization, 
respectively. For long wavelengths, = ^k^, where 
is the bare sound velocity in the direction k. Hi and H2, 
respectively, arc the spin-phonon interactions that arise 
from the first- and second-order terms of Eq. (2), 

^1 = E f^k'^ A, H2 = IY. f^kkUkv4-k' , (5) 

k kk' 

where Ak = 6k + 6^k ^^"-^ 
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Here, M is the mass of the chromium ion (note that 
the same symbol wiU be used for magnetization later; 
the context will clarify the meaning) and N is the to- 
tal number of chromium ions. In addition to nearest 
neighbor exchange, CdCr204 has weak further- neighbor 
exchanges and DM interactions. While the effect of the 
former are expected to be negligible, the DM interactions 
play a crucial role in the low-temperaturc ordered state 
of CdCr204,^° as we will see later. 

Useful simplifications occur in the Hamiltonian given 
by Eq. (3) on noticing that the direct exchange domi- 
nates the magnetic interaction between neighbouring Cr 
spins. Thus, we have Jy = Joe^"^'^ , where a now de- 
termines the strength of the magnetoelastic interactions. 

The effect of the spin-lattice coupling on the spins in 
cubic spinels in similar context was discussed by Tch- 
ernyshyov et al.^ They considered the Hamiltonian given 
by Eq. (1), expanded Jij to linear order {dJij /dTit), 
and showed that the spin degeneracy is lifted through a 
Jahn- Teller-type mechanism selecting a collinearly spin- 
ordered state, and at the same time, distorting the lat- 
tice. On including the second-order term d^Jij/dHfj, 
the potential energy for the bond (ij) is given by 
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where K is the stiffness constant of the bond and J', J" 
are the first- and second-order derivatives of Jij as indi- 
cated in Eq. (2). Integrating out the displacement vari- 
ables (Sij) within the independent-bond approximation, 
we have: 



'2{K + J"S^ ■ Sj 



(8) 



Thus, in the limit K ^ J', J", it is clear from Eq. 
(8) that the second-order term docs not destabilize the 
coUinear spin ordering but renormalizes the different or- 
der parameters at sub-leading order, leaving the general 
structure of the phase diagram intact. (Similar results 
are obtained in a more detailed analysis. "'^^) It is very 
important to note that this is not the case for the renor- 
malization of the phonon spectrum. 

The effect of the spin-phonon interaction on the 
phonon spectrum is effectively studied using perturbation 
theories suited to specific regimes, namely, the paramag- 
netic and the magnetically ordered regimes. We address 
these regimes in the following sections. 



III. THE PARAMAGNETIC PHASE 

In the absence of spin-phonon coupling, in the temper- 
ature regime < T < \Qcw\: the spin correlations tend 
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FIG. 1: (Color online) Change of the sound velocity vs tem- 
perature ai H = for the cl mode in CdCr2 04 for an ul- 
trasound frequency of 107 MHz. The inset shows the region 
around Tn in more detail. Data for 8 T (middle curve) and 16 
T (lower curve) applied along the [111] direction are shown as 
well. For the non zero magnetic fields, the temperature sweeps 
up and down are shown. Notice that the transition temper- 
ature is reduced on increasing the field. All data curves are 
arbitrarily shifted along the y-axis for clarity. 



toward a form that falls off with a fixed integer power of 
distance and a characteristic dipolar angle dependence. 
These correlations lead to a non trivial structure factor 
in the neutron scattering in this cooperative paramag- 
netic phase. ^ Furthermore, in this phase, the spin dy- 
namics has an unusual behavior, ^^'^^ characterized by an 
emergent universal (independent of exchange strength) 
timescale, Tg ^ h/ckBT, c ~ 0(1), that controls the 
long-time dynamics for generic wave vectors. It is known 
that such behavior survives in the presence of weak spin- 
lattice coupling. However, the behavior of the phonons 
has not been studied so far. 

Figure [1] shows the temperature dependence of the 
sound velocity of the cl mode measured below 240 K 
at zero magnetic field. It exhibits a softening below 120 
K followed by a minimum at approximately 13 K and 
a jump like anomaly at Tjv, which is accompanied by a 
small hysteresis signaling the magneto-structural transi- 
tion. The jump shifts to lower temperatures in an ap- 
plied magnetic field [see Fig. 1 (inset) and also Fig. 
3]. This lowering with increasing magnetic field is un- 
derstood from the Clausius-Clapeyron equation for the 
magnetic systems: = Magnetization mea- 

surements (Fig. 7) show AM < as one goes from the 
high-temperature to the low-tcmperature phase. Also, 
AS* < as the system moves into an ordered state from 
the frustrated paramagnet. Hence, < as seen in ex- 
periment. The gradual decrease of the sound velocity oc- 
curs at temperatures corresponding to the paramagnetic 
state of the spins. Note that the infrared reflectivity spec- 
tra exhibit a phonon softening in the same temperature 
range. 

Using the Hamiltonian introduced in Eq. (3) we can 
make quantitative estimates of the phonon self-energy 
correction due to coupling with the spins. In the param- 
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agnetic phase, by comparing the time scales of the acous- 
tic phonons and the spins, we find that the spin dynam- 
ics is fast, and, hence, the spins may be integrated out. 
Here, it is worthwhile to note that the opposite hmit is 
obtained in case of the optical phonons, which is relevant 
for the above-mentioned infrared measurements. ^'^'^^ In- 
tegrating out the spins results in an effective interaction 
among the acoustic phonons given by the effective Hamil- 
tonian. 
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Here (• • • ) denotes thermal averaging over the spins [with 
respect to the unperturbed Heisenberg Hamiltonian Hsp 
given by Eq. (4b)] and ((• • • )) stands for the connected 
correlator: {{AB)) = (AB) - {A){B). Also /3 = 1/KbT, 
where Kb is the Boltzmann constant. Since the timescale 
for the energy exchange between the phonons and the 
spins is inversely proportional to the spin-phonon scatter- 
ing rate, this perturbation calculation is valid only when 
the time scale associated with the energy exchange be- 
tween the spins and the acoustic phonons is the longest 
timescale in the problem. At present, although there is 
no estimate of this timescale, our results provide some 
justification for this assumption a posteriori?^ 

We start our calculation by writing the Matsubara- 
Green's function for the phonons, CJ''(q, t — r') = 
— {T{A^{t)A_^{t')}) . and find the contribution to the 
phonon self-energy to the lowest order in spin-phonon 
coupling. The dressed phonon propagator is given by 



(zao'-K)'-2c^"s(ir!„,q) 



, (11) 



where q is the wave vector, 51„ = 2-Kn/ (3 is the bosonic 
Matsubara frequency and S(if2„,q) is the phonon self- 
energy, which is given by 



E(i^]„,q) = Vq, 



(12) 



where Viqq is given by Eq. (10). From this we find the 
leading-order change of the phonon frequency as 



Awq = ( Wl + ^R-e [S(il^„,q)|.o„^^„+.r,] - 1 



Re[E(ir!„,q)| 



(13) 



Hence, the fractional change in the velocity of sound is 
given by 



lim 



g->0 CJ 
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FIG. 2: (Color online) The best fit (dashed line) of the zero 
magnetic field experimental data (solid line) in the paramag- 
netic regime and determination of a. The non harmonic part 
present in the experimental data (in Fig. 1, for H = 0) has 
been subtracted (using Ref. [29], see text for details). The 
data have been fitted in the regime of T = 20 — 200 K using 
Eq. (16). The variation in a (calculated from these fitting) 
is due to the uncertainty in the value of J (J = 0.80 -1.00 
meV). Note the good quality of the fit down to T <C [Ocwl. 



As shown in Appendix A 1, for the [111] direction of the 
sound, Eq. (14) becomes 
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where C'^p'" and f'^P™ are the unperturbed magnetic spe- 
cific heat and average energy, respectively, of the sin- 
gle spin for a nearest-neighbor Heisenberg AF on a py- 
rochlore lattice [Eq. (4b) with h = 0] . Inclusion of non- 
linear terms into the phonon Hamiltonian (9) (as well 
as probing sound velocity in directions other than [111] 
[see Appendix Al]) may modify the coefficients in front 
of f P'" and C^P™T in the above expression, but we ne- 
glect the corresponding effects. In the cooperative para- 
magnetic phase of the classical Heisenberg AF it was 
shown^® that the energy and specific heat can be cal- 
culated quite accurately using a single tetrahedron ap- 
proximation. Thus, we have [using direct exchange for 
the spin-spin coupling (J = Jqc^"^)] 
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where and C^^^ are the energy and specific heat per 
tetrahedron. To compare with experiments, it is useful 
to factor out JqS^ from the expressions of 6'^°^ — TC^'^^ 
by expressing them as a function T /{JqS"^) = Tg. Thus 



5 



< -10 





■ 1 1 1 ■ 1 ■ 

CdCr^O^ 

H II /f II t; II [111] 
OT 




^^^^.^^ B = 16T 



5 10 15 20 
T(K) 



FIG. 3; (Color online) The sound velocity vs. temperature 
at different fields. Notice that the upturn is suppressed on 
increasing the field (see text for details). The data curves are 
arbitrarily shifi^ed along the y-axis for clarity. 
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FIG. 4: (Color online) Best fit (solid line) of the low- 
temperature data (symbols) (in zero magnetic field) using Eq. 
(18). The data clearly show an ~ behavior as expected 
from the theory of scattering of magnons and phonons. 



Eq. (16) becomes 

V - 18M«)2 J' ^^^^ 

Xet ^Tct 

where E and C are functions of the single parameter 
To, which can be calculated exactly in the absence of a 
magnetic field. ^® The explicit expressions are given in 
Appendix B. The factor 1/2 is multiplied in Eq. (16) 
because each spin is shared by two tetrahedra. 

In the expression given by Eq. (17) the numerical val- 
ues of all parameters, except a are independently known. 
Hence, we can use a single-parameter fit describing the 
experimental data to get an estimate of a. However, 
there is a non-harmonic contribution to the temperature 
dependence of the sound velocity clearly seen in Fig. 1 
at T above 120 K. This contribution is superimposed 
on the sound velocity renormalization arising from the 
spin-phonon coupling. Applying a well-established pro- 
cedure, we have subtracted the non harmonic contribu- 
tion from the experimental velocity data of Fig. 1 us- 
ing an empirical equation from Rcf. 29. The fit of the 
corresponding experimental data is shown in Fig. 2 for 
a = 13.7 ± 0.6 A . (The error bar is mainly due to the 
uncertainty in the determination of the exchange cou- 
pling Jq.) Aguilar ei al?^ and Kant et ai?"^ measured 
the shift of the infrared-active optical phonon. From their 
optical-phonon data, we find a = 10.97±1.24 A . Thus, 
the coupling constants derived from optical and acoustic 
phonons are quite close to each other. 

IV. THE MAGNETO-ELASTIC TRANSITION 

The hysteresis in the sound velocity at T/v confirms 
the first-order type of the phase transition, previously 



also observed in magnetic-susceptibility measurements.^ 
The sound attenuation (not shown) increases on lower- 
ing the temperature and shows a peak like anomaly at 
T]S!?^ The strong anomalies in the acoustic properties at 
the transition reflect the crucial role of the spin-lattice 
coupling in this compound. A particularly interesting 
feature is the smooth upturn (from about T = 13 K) 
in the sound velocity as a precursor to the actual jump 
at the magnetoelastic transition (see Fig. 1). Similar 
anomalous features have also been seen in specific-heat 
measurements.^^ This upturn is suppressed with increas- 
ing magnetic field (Fig. [3]). At present there is no clear 
understanding of its origin. The fact that the upturn 
is suppressed by the magnetic field indicates that it re- 
sults from spin fiuctuations. This suggests the presence 
of near-critical modes in the vicinity of the magnetoelas- 
tic transition that are ultimately cut off at some length 
scale leading to the actual first-order transition. 

A Landau free-energy analysis based on the interaction 
of the magnetic modes and strain fields^^ qualitatively 
accounts for the actual jump of the sound velocity (the 
jump being proportional to the magnitude of the order 
parameter on the low-temperature side). However, such a 
phenomenological theory is highly qualitative due to the 
large number of terms allowed by various symmetries. 



V. THE LOW-TEMPERATURE ORDERED 
PHASE 

In the low-temperature state, the lattice is distorted, 
and the spins are ordered in a long-pitched spiral. Ow- 
ing to the lattice distortion, the erstwhile [111] direction 
is no longer a symmetry axis. Thus, it is difficult to 
extract the magnetoelastic coupling constant (s) in this 
region and to check their consistency with that obtained 
in the paramagnetic regime. 
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However, in this regime, an analysis based on the 
magnon-phonon scattering in a cohinear AF gives in- 
sights into the temperature dependence of the sound ve- 
locity. We can no longer integrate out the spins due to 
the presence of the low-frequency long-wavelength spin 
waves characteristic to the broken spin-rotation sym- 
metry phase. Thus, essentially, we have a problem of 
magnon-phonon scattering, and we wish to calculate the 
renormalization of the phonon energy due to these scat- 
terings. Hence, we start again with the Hamiltonian in 
Eq. (1) and use the Holstcin-PrimakofF (HP) approxi- 
mation for the spin-dependent potentials uj^^ and U^l., 
[Eqs. (6a) and (6b)]. Then, using diagrammatic per- 
turbation theory as before, we calculate the rcnormal- 
izcd phonon propagator [as in Eq. (11)]. The calcula- 
tion of the phonon self energy is a problem of evaluating 
the various magnon-phonon diagrams. ■^^^"^'^ After some 
rather tedious but ultimately straight forward diagram- 
matic calculation [outlined in Appendix A 2] , we obtain 
the phonon self-energy. Then, using Eqs. (13) and (14), 
we compute the fractional change in the sound velocity 
(to the leading order at H = 0) 



Av 



(c + KT^ 



(18) 



The coefficients c and K depend on the form of the lat- 
tice and direction of the sound wave, in general, and arc 
hard to determine theoretically, particularly for a dis- 
torted crystal. However, on general grounds we can ar- 
gue that c < 0. This is essentially the contribution arising 
from in Eq. (6b) with the spins in their classical 

ground state. In this case, it is easy to see that c < 
which indicates that the phonon mode will be softened 
at T = due to the coupling with the spins. Figure 4 
shows such a fit, which compares fairly well with the ex- 
perimental data. [A similar calculation for a ferromagnet 
predicts that ^ cx (c -f KT^).] 



VI. MEASUREMENTS IN NON ZERO 
MAGNETIC FIELDS 

We have already seen that the sound velocity depends 
on short-ranged spin correlations. The magnetic field 
affects the spin correlations, and, hence, the sound ve- 
locity is affected as well. Ultrasound measurements as a 
function of magnetic field (at various temperatures) are 
presented in Fig. 5. Besides some low- field features (be- 
low T^r, which we will discuss later), the acoustic mode 
demonstrates a clear softening with increasing field. This 
general trend may be understood qualitatively as follows. 
In the presence of a magnetic field, the free energy of the 
system is given by 
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(19) 



where IFq is the free energy in the absence of the magnetic 
field and x(e) is the magnetic susceptibility (which is a 
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FIG. 5: (Color online) Change of the sound velocity of the 
cl acoustic mode in CdCr2 04 vs magnetic field measured at 
different temperatures. The ultrasound frequency was 107 
MHz. The arrows indicate the field-sweep directions. The 
experimental geometry is H \\ k \\ u \\ [111]. The curves ob- 
tained at different temperatures are arbitrarily shifted along 
the y-axis for clarity. 



function of the strain field e in the presence of magneto- 
elastic coupling). This gives an additional contribution 



to the elastic constants Ci 
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and this is observed in the experiments (see Fig. 6). 
A more microscopic consideration (for the paramagnetic 
phase) developed in Sec. HI suggests that the fractional 
change in the sound velocity in the presence of a magnetic 
field is given by 
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Here, and f ^ are the scaled [see discussion follow- 
ing Eq. (16) and also Appendix B] exchange parts (see 
below) of the specific heat and magnetic energy, respec- 
tively, i.e., 



NKbT^ 



Sjf))h (22) 



and 



pTot _ ^pspin 



(23) 



where, as before, the factor 2 arises because each tetra- 
hedron contains two spins. However, unlike the case of 
zero magnetic field, these quantities cannot be expressed 
(even within the single tetrahedron approximation) in 
closed analytic form. Hence, we resort to Monte Carlo 
simulations of a Heisenberg magnet on a single tetrahe- 
dron to calculate these expressions. We use a heat-bath 
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FIG. 6: (Color Online) A fit (solid line) of the velocity varia- 
tion (symbols) with magnetic field from experiments at T = 8 
K. 

Monte Carlo algorithm and find that we can fit the sim- 
ulation data in the following phenomenological form: 

(£r-ToCr)=<^(To)(^)\ (24) 

where <f>(To) is a function of only the scaled temperature 
introduced before [To = T/{JoS'^)] and H is the mag- 
netic field measured in Tesla. Using this phenomenolog- 
ical form of the function [Eq. (24)] in Eq. (21), we get 

This is the expression for the dependence of the sound 
velocity on the magnetic field in the cooperative para- 
magnetic regime. This expression is in accordance with 
our earlier expectations that Av/v (x H^. While a deter- 
mined from these data is broadly consistent with our ear- 
lier estimation, this agreement is no longer quantitative- 
the data analyzed here are at temperatures too low for 
our perturbation theory. 



A. Low-temperature ordered phase and role of the 
DM interaction 

For temperatures below TjVi there is a kink like 
anomaly at about 4 T accompanied by hysteresis, see 
Fig. 5. Magnetization measurements in this temperature 
and magnetic-field region (Fig. 7) also exhibit a similar 
hystcrctic behavior in the temperature range below Tn 
with an anomaly at a magnetic field of about 4 T. This 
is clearly visible for the M/H vs H data shown in the 
inset of Fig. 7. Above Tat, the magnetization exhibits a 
linear magnetic-field dependence. 

Chern et al}'^ showed that the long-ranged spiral 
results from weak DM interactions. Their presence 
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H(T) 

FIG. 7: (Color online) Magnetization measured for magnetic 
fields applied along the [111] direction at temperatures of 2, 4, 
and 8 K. Zero-field cooled curves are shown. The inset shows 
the same data presented as M/H vs H. The arrows indicate 
the field-sweep direction. 

[Hdm = Z]{jj> Dij • Si X Sj, with non-collinear Dy) 
breaks the global spin-rotation symmetry to discrete lat- 
tice symmetries and results in an energy gap (per tetra- 
hedron) of the order of A « 0.7 K (obtained from 
ab initio calculations^^). The Zeeman term, on the 
other hand, favours a four-sublattice canted spin struc- 
ture. A first-order transition between the spiral and the 
canted states"'^^'"'^^ is responsible for the observed low- 
temperature anomalies. An estimate of the magnitude 
of the magnetic field at which this transition occurs may 
be obtained by comparing the energies of the spiral and 
the canted states. The Zeeman-energy gain (per tetra- 
hedron) in the canted state is given by Et = aH^ with 
a = -0.0431 {K/T'^) and H in Tesla.^^ Thus, a rough 
estimate of the magnitude of the magnetic field at which 
the transition occurs is ~ 4.3 T. (We have assumed 
the case of a strongly first-order transition where the 
undistorted spiral changes suddenly to the canted state 
at H = He- The hysteresis in the sound- velocity and 
magnetization data justifies our assumption.) 

Below He, in the spiral state, the spin correlations do 
not change appreciably with the magnetic field leading 
to an almost constant sound velocity, as seen in exper- 
imcnt(Fig. 5). Above (He), in the canted state, the 
spin correlations, however, are sensitive to the magnetic 
field, and the sound velocity changes with magnetic field 
{Av/v (X H^). Our estimate of He matches fairly well 
with the value of the field at which anomalies are seen in 
the experiments (Figs. 5 and 7). But it is not clear why 
the characteristic field of 4 T is somewhat lower than 5.7 
T, where similar anomalies have been reported in ESR 
measurements . ^ ^ 



B. High magnetic field regime 

We have also performed pulsed magnetic-field experi- 
ments in CdCr204 in order to study the acoustic behav- 
ior in the region of the magnetization plateau. Results 




20 40 60 

H(T) 

FIG. 8: (Color online) Change of the sound velocity for the cl 
mode in CdCr2 04 measured in pulsed magnetic fields at dif- 
ferent temperatures and an ultrasound frequency of 81 MHz. 
The experimental geometry is H \\ k \\ u \\ [111]. 



transition. It is worth noting that the magnetization 
only exhibits a smooth kink like anomaly at this phase 
transition. A transverse spin order, which is equivalent 
to a Bose-Einstein condensation of magnons, is predicted 
from quantum-fluctuation theory at magnetic fields just 
above the plateau. 

The first-order phase transition at about 28 T could 
be resolved by magnetization measurements up to tem- 
peratures slightly above Tm-^ Our pulsed- field ultrasound 
measurements performed at 7.8 K (Fig. 8) also clearly re- 
veal an anomaly in the sound velocity at this field. The 
anomaly at 28 T is somewhat smoother than at lower 
temperatures but still clearly evident. The hysteresis sur- 
vives, and the total change in the sound velocity is even 
larger at higher temperatures approaching 2 % between 
and 60 T. Further temperature increase suppresses the 
anomaly at 28 T leading to a broad minimum in the 
sound velocity at about 50 T (see the lower curve in Fig. 
8). Note that the hysteresis is still observable at 12.5 K. 



for the sound velocity in magnetic fields up to 63 T at 
temperatures of 1.4, 4.7, 7.8, and 12.5 K are shown in 
Fig. 8. At 1.4 and 4.7 K, the sound velocity decreases 
first with the characteristic behavior, then demon- 
strates a jump at the magnetic field where the mag- 
netization plateau appears. This anomaly corresponds 
to a first-order phase transition from the four-sublatticc 
canted spin structure to a collinear spin configuration 
with three spins up and one spin down at each Cr^"*" tetra- 
hedra. Previously, a large magnetostriction has been re- 
ported at this phase transition.^ A cubic crystallographic 
structure has been suggested from high-field x-ray ex- 
periments at the plateau state. '^^ Recent elastic-neutron- 
scattcring experiments in pulsed magnetic fields'" showed 
that the magnetic structure at the half-magnctization 
plateau phase has a cubic P4332 symmetry. It has been 
proposed^^'^^ that the lattice distortion stabilizes the 
three-up one-down collinear spin configuration. Once 
again, the lattice distortion complicates the theoretical 
analysis of this regime. Indeed, there is only a slight 
change of the sound velocity within the magnetization- 
plateau range. This is because the spin correlations arc 
locked at fixed values within this plateau. The plateau 
terminates at approximately 58 T, confirmed by a sharp 
anomaly, i.e., an abrupt decrease in the sound veloc- 
ity. This anomaly corresponds to a phase transition to 
a non-colhnear canted spin configuration. ^'^'■^^ It is inter- 
esting to note that the hysteresis in the sound velocity 
not only takes place at the first-order phase transition 
around 28 T, but also spreads along the whole plateau 
range up to 58 T, showing a complicated interplay be- 
tween the spin and the lattice degrees of freedom within 
the magnetization plateau. The sound-velocity change, 
which takes place between 58 and 63 T, is even larger 
than the anomaly at 28 T. The highest applied mag- 
netic field of 63 T is not sufficient to detect the complete 
sound-velocity change. No hysteresis has been detected 
at 58 T, pointing to a second-order type of this phase 



VII. CONCLUSION 

In this paper, we have successfully characterized the 
strength of the magneto-elastic coupling in CdCr2 04 and 
have been able to explain the general features except at 
highest fields. Wc would like to mention that, in the 
current paper, we deal with the acoustic cl mode where 
various deformations are involved (see the definition of 
Cl mode above). This fact complicates the symmetry 
analysis of the obtained data. In addition, below T^r, in 
the tetragonal phase, there are three types of domains 
corresponding to an elongation of the c axis, and all of 
them contribute to the acoustic cl mode. In this re- 
gard analysis of the sound mode along one of the axes 
of the cubic/tetragonal crystal may yield more theoreti- 
cally tractable results. A detailed analysis of these other 
modes and their relation to the different elastic constants 
are given in Ref. 17. 

We also note that hydrodynamic calculations suggest 
that, in the regime |6cw| ^ T ^ Tm (where Tm is 
the melting temperature of the crystal), Av/v — )■ K, 
where TiT 7^ is a constant. Thus, in this limit, there 
is a constant shift in the sound velocity compared to an 
iso-structural compound without spin-phonon coupling. 
This result is rather interesting and requires further in- 
vestigation. 

In general, the sound waves, being hydrodynamic 
modes, are very robust, and they couple to a large num- 
ber of low-lying excitations. It is relatively straightfor- 
ward to measure the sound velocity with high enough 
accuracy, and this is especially true for the present case 
of magnetoelastic transitions. However, since they them- 
selves do not go critical and also because they generally 
interact with most low-energy modes (mentioned above) 
disentangling microscopic information from them is less 
straightforward. 

To summarize, we have presented a magneto-acoustic 
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study of the frustrated spin system CdCr204. Strong 
sound- velocity anomalies have been observed at the mag- 
netic phase transitions in CdCr204. We have been able 
to characterize the spin-strain coupling, which is crucial 
and determines the underlying physics of this compound, 
reducing the geometric frustration effect and lifting the 
degeneracy in the system. The detailed comparison of 
experiment and theory in different regimes suggests that 
the dominant part of the variation of the sound velocity 
at low temperatures is due to the spin-lattice interaction, 
which in turn, may be modeled as an exchange-striction 
phenomenon. 
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Appendix A 

In this appendix, we give an outline of the perturbative 
calculation for both the high-temperature paramagnetic 
phase as well as the low-temperature ordered phase. 

1. High-temperature paramagnetic phase 

Here, we derive Eq. (15) from Eq. (14) for the sound 
velocity in the [111] direction. From Eqs. (10) and (12), 
we find that the phonon self-energy is given by 

S(za.q) = V^^ = {U^^) - /3((C/(i)c/i^^)), (Al) 
where, from Eqs. (6a) and (6b), we have 



2MNuj' 



E E«(S. • S,) (S, • S^))){e^^-' - e^^-n 



{ij) {Irn) 



(ij) 



(A2) 
(A3) 



Taking the limit given by Eq. (14) we have (for longitudinally polarized phonons) 



Av 



V 2MN{viy 



Jo dS^ 



(ij) 



, (A4) 



where, 5 is the distance between two neighbouring Cr spins, which is same for all directions (%). In a paramagnetic 
state and a general direction of the sound velocity (q), we have 



Av 



V 2M(uy 



0\2 



^(q)f ^P'- - fiS ) S(q)C^P'"r 



Jo dS^ 



Jo 86 



(A5) 



r 



In the last expression, we have used 



N 
1 

N 



Jo Si • Sj 

iv) 



(A6) 
(A7) 



For q [111], we have (from the fact that S ■ q = 
exactly half of the bonds), 

^(q) = B{q) = I 



t for 



(A8) 



Hence, we get Eq. (15). 

2. Low-temperature magnetically ordered phase 

Here, we outline the calculation for the fractional 
change of sound velocity in the low-temperature ordered 
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state due to magnon-phonon scattering. We assume a 
collinear two sub-lattice Neel order. While the differ- 
ences in the lattice structure may lead to variations in 
different non-universal pre factors, the temperature de- 
pendence does not change as long as the magnon disper- 
sion remains linear. Similar calculations (not shown) can 
be performed for the ferromagnetic case as well. We start 
by defining two species of HP boson operators, 

Sublattice 1 : 5^ = V2Sa, = ^/2Sa\ ^ S - flta, 

Sublattice 2 : = V2Sc, = ^/25c^ ^ S - c)c. 

(A9) 

The Hamiltonian is given by 

H = Hp + Hsp + H1+H2, (AlO) 

where 

Hp = ri^c^06[5k, (All) 

k 

Hs, = -^^!^-E-£(4"k + /3i/?k).(A12) 

k 

(z is the coordination number; uj^ — Csk is the AF 
magnon frequency.) a and /3 are the Bogoliubov rotated 
bosonic operators given by 



Ok 




cosh 9k sinh 6'k 




Qk 






sinh 6'k cosh 0k 







(A13) 



and 

tanh26lk = -^^^ . (A14) 



Also 

H, = E^k*Uq-Xk,q-*q, (A15) 
k,q 

k 

+ ^ E ^kA_k'*Lk' + q-A/k,k',q-*q,(A16) 
k,k',q 

where ^'j^ = [aj^, /3_k], ^k = ^k + &lk i'' ^^'^ phonon 
displacement operator and 

(A17) 

and A^k.q and A/k.k'.q are 2x2 matrices 

Xk,q = Ck+q • Ak.q • Cq, (A18) 
A/ick'.q = Ck-k'+q ■ C'k,k',q ■ Cq, (A19) 

where, 

/cosh 0k sinh 0k \ (A20) 
^ \ smh0k cosh 0k / 

and 



Ak,q = - 



JoSa 
a^JoSn 



2MNuj° 

^ s 



E (1- 



{S ■ ek) 



^k,k' q — / rr 



k' S 



E (1- 



■ik-i5\ 



1 - e 



— ik'-i5 



1 6*1'' 
gl(k+q)-<5 2 

{6 ■ ek)(i5 • e_k') 



1 

„j(k-k'+q)-5 



(A21) 
(A22) 



The phonon Green's function is given by Eq. (11) and 
the self-energy consists of two parts 

I](q,if]„) = I]i(q,il]„) + E2(q,«f^«). (A23) 

Defining the bare matrix Green's function for the 
magnons as Gj^.(k,T) = -(r{*U^)V'( V)}> or its 



Fourier transform, 

G™(k,*aO - ( "1 ) , (A24) 

we find that the contributions to the self-energy from Hi 
and H2 are given by 
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^2{^n„, q) = r(°) + 1 ^ (AA^^q,^, + AA^'q.^q.qO Gj;(zO-, q'). (A26) 
q' 



Now, the frequency and the fractional change of sound 
velocity are given by Eqs. (13) and (14). In general, the 
integrals in Eqs. (A25) and (A26) are hard to calculate 
for a general direction of q. While this can be done nu- 
merically, here, we only concentrate on their dependence 
on temperature, which gives the result given by Eq. (18). 



I 

neighbour classical Heiscnberg AF (with spins of unit 
magnitude) on the corners of a tetrahedron. The tem- 
perature is scaled by l/(Jo5^), i.e., Tq = T/JqS^, to 
compare with experiments. 



Appendix B 



^Tct 



Here we give the explicit expressions for £ and 
C , which are the energy and specific heat of nearest- 



J 



Xct Tct 

£ and C in Eq. (17) can be obtained from Ref. 28 and are given by 



-^■Tot 



£ = 

:Tet 3 



c 



e2 




Tot , 

Z 732^ 


2/To 


if''- 


\/32^ 





- To 3 - 4e-2/^» + +61- e^^/To 



-Tet 



where Z""" is the scaled partition function of the tetrahedron, which is given by^® 



2Erf 



(72771;) - Erf [^8j%) - V ^ (1 - e"'^''") (3 - 2e-^/^« + e"^/^") 



Here, Erf (cc) is the error function defined as Erf{x) — 1/4/71 Ji^ e * dt. 

r"^ 



(Bl) 
(B2) 



(B3) 
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